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Abstract 

Hydrodynamic equations for ideal incompressible fluid are written in terms 
of generalized stream function. Two-dimensional version of these equations is 
transformed to the form of one dynamic equation for the stream function. This 
equation contains arbitrary function which is determined by inflow conditions 
given on the boundary. To determine unique solution, velocity and vorticity 
(but not only velocity itself) must be given on the boundary. This unexpected 
circumstance may be interpreted in the sense that the fluid has more degrees of 
freedom, than it was believed. Besides, the vorticity is less observable quantity 
as compared with the velocity. It is shown that the Clebsch potentials are 
used essentially at the description of vortical flow 
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1 Introduction 



In this paper we write hydro dynamic equations for ideal fluid in terms of Clebsch po- 
tentials PQEj- This representation admits one to integrate hydrodynamic equations 
and to obtain effective description of rotational stream flow. We have discovered 
unexpected fact that the effective description of the rotational flow cannot be car- 
ried out without introduction of the Clebsch potentials. In general, the problem of 
irrotational stream flow and that of rotational stream flow are stated differently. Let 
us show this in example of stationary two-dimensional flow of incompressible fluid. 
In this case one of Clebsch potentials may coincide with the stream function if). The 
irrotational flow is described by the equation for stream function if) 

where indices mean corresponding partial derivatives. For the rotational stationary 
two-dimensional flow the equation transforms to the form 

^+^ w = nW (i-2) 

where Q (if)) is some function of argument if), describing the character and intensity 
of vorticity. The obtained equation distinguishes from (jl.lj) in additional term. 
Both equations (jl.lj) and (jl.2j) are the elliptic type equations. To obtain a unique 
solution of equation (jl.lj) . it is sufficient to give value of the stream function if) 
on the boundary X of the volume V, where the flow is considered. The same is 
valid for equation (jl.2j) . provided the form of the function Q (if>) is known. But the 
value of if) on £ does not admit one to determine the form of the function Q (if>). 
For determination of Q (if)) one needs some additional information, given on the 
boundary E. 

To state the problem of the rotational stream flow, it is necessary to determine 
what kind of information on the boundary is sufficient for derivation of the unique 
solution of equation (jl.2j) . Theory of rotational flows cannot exist without statement 
of this problem. Unfortunately, we have not seen in literature a statement of the 
problem of the rotational stream flow, we know only single exact solutions |Hj- It 
means that the theory of rotational flows does not exist. 

Equation (jl.lj) is a special case of equation (jl.2j) . when Q(if>) = 0. Equation 
(jl.2j) is obtained from the conventional hydrodynamic equations for two-dimensional 
stationary flow of incompressible fluid, which have the form 

U X +Vy = 0, dy (UU X + VUy) = 8 X (UV X + VVy) (1.3) 

where u and v are velocity components along x-axis and y-axis respectively. Intro- 
ducing the stream function if> 

u = -ifj yJ v = ip x (1.4) 

we satisfy the first equation (jl.3j) identically. The second equation (jl.3j) takes the 
form 

d(u,if>) 

j^ = 0, *^„ + ^ w (1.5) 
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Relation ()1.2|) is the general solution uj = of ()1.5|) . where Q is an arbitrary 

function of ip. 

The goal of the present paper is the statement of the problem of the rotational 
stream flow. The problem is solved effectively in terms of the generalized stream 
function (GSF) which has several components. GSF is a generalization of the usual 
stream function ip for two-dimensional flow onto a more complicated n- dimensional 
case. Unfortunately, in other more complicated cases the statement of the problem 
of the rotational flow needs a special well developed technique. 

This technique is based on application of Jacobian technique to the descriptions 
of hydrodynamic equations written in terms of Clebsch potentials. We shall refer to 
this technique as the generalized stream function technique (GSF-technique). Use 
of Jacobian technique in application to Clebsch potentials goes back to papers by 
Clebsch |TJ |2] , who obtained his potentials dealing with Jacobians. In contemporary 
hydrodynamics the Clebsch potentials are considered formally without connection 
with the Jacobian technique 0]. Clebsch potentials have also another names (La- 
grangian coordinates, Lagrangian variables [S], labels of fluid particles). There are 
several version of representation of Clebsch potentials. Elements of Jacobian tech- 
nique was used by many authors ( El El IS] and many others) . But in solution of 
the problem of streamline flow the Jacobian technique is not used practically, as well 
as Clebsch potentials. We are interested only in that version, which is connected 
with a use of Jacobian technique (GSF-technique). We use space-time symmet- 
ric version of the Jacobian technique which appears to be simple and effective. It 
seems that the progress in the investigation of vortical flows is connected mainly 
with the developed Jacobian technique used together with Clebsch potentials (GSF- 
technique). Presentation of the GSF-technique can be found in ^U|, as well as in 
[TTj . where it is used for description of the fluid flow in terms of the wave functions. 

Why are Clebsch potentials important in the statement of the problem on vortical 
stream flow? The term Q (-0) in rhs of (jl.2J) describes vorticity. According to the 
Kelvin's theorem on circulation the vorticity is "frozen in the fluid", and vorticity 
travels together with the fluid. Clebsch potentials (labels of the fluid particles), as 
well as the generalized stream function (GSF) are also "frozen in the fluid". They 
also travel together with the fluid and vorticity. As a result the Clebsch potentials 
(labels) and the generalized stream function (GSF) appear to be an effective tool of 
the vorticity description. 

A rotational fluid flow has some properties which are absent in irrotational flows. 

1. Consideration of the fluid displacement and a use of Clebsch potentials describ- 
ing this displacement is essential in rotational flows, because this displacement 
transfers vorticity, which influences the velocity field, whereas such a transport 
is of no importance in irrotational flows, where the vorticity vanishes. 

2. The boundary conditions for a rotational flow contain more information, than 
the boundary conditions for irrotational flow. This additional information is 
information on vorticity, which conserves in any ideal fluid. Although the field 
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of vorticity 



oj = V x v 



is determined by the velocity field v, but it is valid only inside some 3- volume 
V, where vorticity can be determined as a result of differentiation of the veloc- 
ity field. On the boundary X of 3- volume V one can calculate only component 
of uj along the normal to S. Components of vorticity u) tangent to S must 
be given additionally. It means that additional (as compared with the irrota- 
tional flow) degrees of freedom appear in the rotational flow, and additional 
information in boundary conditions is necessary for their description. 

The first property is analyzed in jlUj . We present it here also. It is common 
practice to think that the problem of streamline flow can be solved by consideration 
of only Euler system of equations 



dp 
dt 



V(pv) = (1.6) 



<9v , _ N 1_ 2 dE 

_ + (vV)v = — Vp, P = Pjt L7 

at p op 

— + (vV)S = (1.8) 

where p is the pressure, and E = E(p, S) is the internal energy of an unit mass 
considered to be a function of the mass density p and the entropy S. The internal 
energy E = E(p, S) is an unique characteristic of the ideal fluid. Displacement of 
the fluid particles (i.e. their trajectories and the law of motion along them) in the 
given velocity field v is described by the equations 

^ = v(t,x), x = x(t,£) (1.9) 

where v(i, x) is a solution of the system (jl.fij) - (jl.8|) . The equations (jl.fi j) - (jl.9j) form 
the complete system of hydrodynamic equations, but the form of this system is not 
quite consistent, because dynamic variables p, v, S in (jl.6j) - (jl.8j) are functions of 
t, x, whereas x in (jl.9j) is a function of t and of the label £ of the fluid particle. 

Let £ = £(£, x) = {£ a (i, x)}, a = 1,2,3 be three independent integrals of equa- 
tions (jl.9j) . Values £ of three integrals can label fluid particles (Lagrangian coordi- 
nates). Then labels £ satisfy three equations 

§ + (vV)£ = 0, (1.10) 

which are equivalent to three equations (jl.9|) . The system of equations (jl.6j) - 
(jl.9j) is hardly perceived as a whole dynamic system, because the system (jl.6Jl - 
(jl.8j) is closed and its dynamic variables p, v, S are functions of t, x, whereas the 
dynamic variables of (jl.9j) are functions of variables t, But equivalent system 
(jl.(jj) - (jl.8j) . (jl.lUj) of equations for variables p, v, S, £, which depend on t, x, is the 
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whole system of dynamic equations. This system is obtained from the variational 
principle, whereas the Euler system (jl.6|) - (jl.8j) of hydrodynamic equations can be 
obtained from the variational principle only for the case of irrotational flow. The 
system (ll.fij) - (|1.8j) is a closed subsystem of the whole system of dynamic equations 
ffTTfijl -fTK ]) . ifTTTTTjl . On the foundation of closure the Euler system (fTfijl - fTSjl is 
considered conventionally as the complete (full) system of hydrodynamic equations, 
whereas in reality the Euler system (jl.6|) - (jl.8|) is only a curtailed system, i.e. only a 
part of the complete (full) system of hydrodynamic equations formed by equations 
(ll.fi|) - ()1.8|) . (jl.K)j) . If we work with the Euler system only, we cannot integrate 
it, in general. If nevertheless we integrate it in some special cases, the arbitrary 
functions of Clebsch potentials £ appear in integrated dynamic equations. If we 
use the complete system (|1.6|) - (jl.8j) . ()1 . lUj) . we can always integrate it and reduce 
the number of dynamic equations. This integration is accompanied by appearance 
of three arbitrary functions g (£) = {g a (£)} , a = 1, 2, 3 of Clebsch potentials £ in 
dynamic equations. These functions contain full information on initial and boundary 
conditions for the fluid flow. 

The integrated dynamic equation have different form for the irrotational and 
rotational flows. In particular, if V x g = 0, what corresponds to irrotational 
flow of barotropic fluid, the equations (jl.lO)) . known as Lin constraints ^2j, are not 
dynamic equations in the integrated system of dynamic equations, and the variables 
£ have no relation to calculation of the fluid flow. In general case of rotational flow, 
when V x g / identically, the equations ()1.10|) are dynamic equations for the 
fluid. We show this, as soon as the corresponding mathematical technique will be 
developed. 

In the second section Jacobian technique is considered. The third section is de- 
voted to derivation of hydrodynamic equations of ideal fluid from the variational 
principle. The variational principle for incompressible fluid is considered in the 
fourth section. In the fifth section two-dimensional flow of incompressible fluid is 
described in terms of the generalized stream function. In the sixth section the inflow 
conditions for the stationary two-dimensional flow of incompressible fluid are intro- 
duced. The seventh section is devoted to stationary flow around an obstacle. Some 
examples of two-dimensional stationary flow of incompressible fluid are considered 
in the eighth section. 



2 Jacobian technique 

Let us consider such a space-time symmetric mathematical object as the Jacobian 



J 



d(£ ,£i,£ 2 ,£ 3 ) 
d(x°, x 1 , x 2 , x 3 ) 



det 



i.k 



i,k 



dkii 



dx k 



z,fc = 0,l,2,3 (2.1) 



Here £ = {£ ,£} = {£ , £ 1; £ 2 > £3} are f° ur scalar considered to be functions of 
x = {x°,x}, £ = £(x). The functions {£cb £15 £2; £3} are supposed to be independent 
in the sense that J ^ 0. It is useful to consider the Jacobian J as 4-linear function 
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of variables Cik = dk£,i, i, k = 0, 1, 2, 3. Then one can introduce derivatives of J with 
respect to fc . The derivative dJ/d^ k appears as a result of a replacement of £ 4 by 
x fc in the relation (12.11) . 

i, = 0, 1,2,3 (2.2) 



8 J _ d(Z ,...£ i _ 1 ,x k , £ i+1 ,...£ 3 



For instance 

«9J _ d(**,£i,£ 2> £3) i = 0jl>2> 3 (2.3) 



9£o,i 8(x°, x 1 , x 2 , x 3 ) 
This rule is valid for higher derivatives of J also. 



9€i,k d Za,i 8(x°,x l ,x 2 ,x 3 ) 
_ d(x k ,x l ) djto,^,^,^) _ 8x k 8x l dx k 8x l _ ni 

It follows from (J2HJ) , Q that 



dx k _ d(t ,...£ i _ 1 ,x k ,£ i+1 ,...Q 

a^o,.-^-!,^,^!,--^)^ ,^ 1 ,^ 2 ,^ 3 ) _ 1 <9J 



(2.5) 



i,Jfe = 0,1,2,3 



and (J2.4j) may be written in the form 

d 2 J 1, dJ 8J 8J 8J 



i, k, I, s = 0,1,2,3 (2.6) 



The derivative dJ/d^ k is a cofactor to the element £ f fc of the determinant (|2.1I) . 
Then one has the following identities 

8 7 81 

^ =8tJ, ^— = SU, l,s = 0,l,2,3 (2.7) 

8 J 8 2 J 

Here and in what follows a summation over two repeated indices is produced (0-3) for 
Latin indices. The identity (j2.8j) can be considered as a corollary of the identity (|2.fijl 
and a symmetry of 8kdi^ s with respect to permutation of indices k, I. Convolution 
of ()2.6|) with 8k, or d\ vanishes also. 

8 2 T 8 s J 

dk 8T7)T = 8f 8f 8f dkd ^m = 0, i,l,s = 0,1,2,3. (2.9) 



Relations (|2.1|) - (J2.6)) are written for four independent variables x, but they are 
valid in an evident way for arbitrary number of n + 1 variables x = {x°, x 1 , . . . x n } 
and£ = {£ ,£}, £ = {£i,£ 2 ,---U- 

Application of the Jacobian J to hydrodynamics is founded on the property, 
which can be formulated as the property of the generalized stream function (GSF) 
£ ={£i, £21 •••> £n} i n t ne (n+l)-dimensional space of coordinates x = {x°, x 1 , x n }. 

On the basis of the generalized stream function £ one can construct n-vector j % 

dJ 

f = m ^r~ ' ? = & P y }i m = const ( 2 - 10 ) 

in such a way that j l satisfies the continuity equation 

Qtf* = (2.11) 

identically for any choice of variables Besides the variables £ are constant along 
any line £ tangent to n- vector j\ and can label this line, because the set of quantities 
£ is different for different lines C. 

In the case of dynamical system (jl.6|) - (JHHJ), (|1.1U|) we have j l = {p,pv}. It 
means that equations (jl.6|) . (|1.10|) are satisfied at any choice of the generalized 
stream function (Clebsch potentials) Substituting the flux vector f = {p, pv}, 
expressed via in remaining hydrodynamic equations, we obtain dynamic equations 
for determination of generalized stream function This procedure is insignificant 
in the case of irrotational flow, when Lin constraints (jl.lOj) are of no importance. 
But it is a very effective procedure in the case of rotational flow, because it reduces 
the number of dynamic equations to be solved. 

The continuity equation 1)2.11)1 is used without approximation in all hydrody- 
namic models, and the change of variables {p, pv} <-> £ described by (12.10)) appears 
sometimes to be useful. In particular, in the case of two-dimensional established 
flow of incompressible fluid the generalized stream function £ reduces to one vari- 
able £ 2 = known as the stream function. In this case there are only two essential 
independent variables x° = x, x l = y, and the relations ()2.1L)|) . ()2.11)) reduce to 
relations 

1 dib , dib du dv , 



Defining the stream line as a line tangent to the flux j 

dx dy 

Jx Jy 



(2.13) 



we deduce that the stream function is constant along the stream line, because ac- 
cording to two first equations (j2.12j) . ip = tp(x,y) is an integral of the equation 
(EH- 
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In the general case, when the space dimension is n and x = {x°, x 1 , . . . x n }, 
£ = {£o> £}> £ = {£ i, &2> • • • f nh t]i e quantities £ = {£ a }, a = 1, 2, ... n are constant 
along the line £ : x = x{t) tangent to the flux vector j = {j 1 }, i — 0, 1, . . . n 

C: — =f(x), z = 0,l,...n (2.14) 

where r is a parameter along the line £. This statement is formulated mathemati- 
cally in the form 

= ^« = m |^- = °> a = 1, 2, ... n (2.15) 

The last equality follows from the first identity ()2.7|) taken for s — 0, I — 1,2, . . .n 
Interpretation of the line (|2.14|) tangent to the flux is different for different cases. 
If x = {x°, x 1 , . . . x n } contains only spatial coordinates, the line (|2.14j) is a line in the 
usual space. It is regarded as a stream line, and £ can be interpreted as quantities 
which are constant along the stream line ( 1.6. cLS cL generalized stream function). 
If x° is the time coordinate, the equation ()2.14j) describes a line in the space-time. 
This line (known as a world line of a fluid particle) determines a motion of the fluid 
particle. Variables £ = {£i,£ 2 ; ■ ■ -£n} which are constant along the world line are 
different, generally, for different particles. If £ a , a — 1, 2, ... n are independent, they 
may be used for the fluid particle labeling. When one of coordinates x is timelike 
the set of variables £ is not perceived as a generalization of the stream function 
if). Nevertheless, we shall use the term "generalized stream function" in all cases, 
because from mathematical viewpoint it is of no importance whether coordinate x° 
is timelike or spacelike. 

Thus, although interpretation of the relation (|2.1(J|) considered as a change of 
dependent variables j by £ may be different, from the mathematical viewpoint this 
transformation means a replacement of the continuity equation by some equations 
for the generalized stream function Difference of the interpretation is of no 
importance in this context. 
Note that the expressions 

f = ™Pb(*)#- = ™Po(£)|^%%% * = °> !> 2 ' 3 > ( 2 - 16 ) 
ot, i a(x u , x l , x 2 , x 6 ) 

can be also considered as four-flux satisfying the continuity equation (j2.11j) . Here 
m is a constant and p (£) * s an arbitrary function of £. It follows from the identity 

f t \ 9(x i , £ lf £ 2 , £3) d(x l , £ 2 ; £3) ~c f (d C C \AC> /n 17\ 

m Po(£) d(r O r l T 2 -3N = m fl(r 'I T 2 T 3 V ^ = / Pottl't* W^l- ( 2 - 17 ) 

The Jacobian technique is very useful for manipulation with hydrodynamic equa- 
tions ()1.6|) - ([1.8|h ([1.10|h For instance, one can integrate the complete system of 5+3 



8 



dynamic equations (jl.6j) - (|1.8|) . (|1.10|) reducing it to the system of five dynamic equa- 
tions written in the form JU] 



S(t,x) = S (£) (2.18) 

P^)=Po (*) |%%# = A>(0 |® (2-19) 

v(t, x) = u(>, 1, 77, 5) = V<p + <f(£) V£ Q - (2.20) 

where So(£), Po(£) an d = a = 1, 2, 3 are arbitrary integration functions 

of argument The quantities tp, r\ are new dependent variables satisfying dynamic 
equations 

^ + ufa & 77, 5) - i[ufa tv, S)} 2 + ^ = (2.21) 

^+ufa€,i7,5)V»7 = -^. (2.22) 

If five dependent variables </?, £, r? satisfy the system of equations Ijl.lOj) . (|2.21|) . 
(12.22)1 . five dynamic variables S, p, v (|2.18j) - (|2.2()jl satisfy dynamic equations (jl.fijl - 
(11.8)1 . Indefinite functions <So(£), p (£), and <? a (£) can be determined from initial 
and boundary conditions in such a way that the initial and boundary conditions for 
variables ip, 77 were universal in the sense that they do not depend on the fluid 
flow [TU]. Further dynamic equations ([l.lOp . (|2.21|) . ()2.22)1 will be derived directly 
from the variational principle. 



3 Variational principle 



The ideal (nondissipative) fluid is a continuous dynamical system whose dynamic 
equations can be derived from the variational principle with the action functional 

m , dx , 



A[x] = /{y(f ) 2 -V} Po (Z)dtdt, (3.i; 



where the fluid particle coordinates x = {x a (t, £)}, a = 1,2,3 are dependent vari- 
ables considered to be functions of time t and of labels (Lagrangian coordinates) 
£ = £ 2 , £3}- dx/dt is a derivative of x with respect to t taken at fixed 

dt ~ d(t,^,^Q ~ d{t,t) a - L > 2 > 6 W 

Po(£) is some non-negative weight function, and V is a potential of a self-consistent 
force field which depends on £, x and derivatives of x with respect to £. m =const 
is a mass of the fluid particle. The value of mass is unessential, and without loss of 
generality one may set m = 1. For the ideal fluid the potential V is such a function 
of £ and dx/ d£ which can be represented in the form 

V(£,dx/dt)=mE(p,S), p = p Q ^)( ^^ \ , S = S (£) (3.3) 
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where the entropy per unit mass S = Sq(£) is a given function of £ determined by 
the initial conditions, p is the fluid density, and E(p, S) is the internal energy of the 
fluid per unit mass of the fluid. The quantity E(p, S) is considered to be a given 
function of its arguments. 

The variational principle (J3.1)) generates the sixth order system of dynamic equa- 
tions for six dependent variables x, dx./dt considered to be functions of independent 
variables t, £. Such a way of description is known as the Lagrangian description of 
the fluid. If the variables t, x are considered to be independent variables, and six 
variables £, d£/dt are considered to be dependent variables, the sixth order system 
of dynamic equations for £, d^/dt arises. It is known as the Eulerian description of 
the fluid. 

The partial choice of labeling of the fluid particles by variables £ is not material 
from a physical point of view. This circumstance is displayed in an existence of the 
relabeling transformation group 

£a-la = ta(£)> D = det II C/^ U 0, «, /? = 1,2,3 (3.4) 

The action (J3.1|) appears to be invariant with respect to the relabeling group, pro- 
vided V has the form (J3.3)) . where S is a scalar and p(£) is a scalar density, i.e. 
under the transformations ()3.4|) p (£) transforms as follows 

Po(0 - Ml), P(0 - Hi) = PiOD (3.5) 
The relabeling group is used in hydrodynamics comparatively recently [T3*l [7J IT4*1 

ma maim ei. 

The relabeling group is a symmetry group of the dynamical system. It may be 
used to simplify a description of the fluid. There are at least two different ways of 
usage of the relabeling group. The first (conventional) way of such a simplification 
is a use of five relabeling-invariant variables v = dx/dt, p, S considered to be 
functions of independent variables t, x. This way of description will be referred to 
as the relabeling invariant description (RID). The variables describing labeling, 
are eliminated. At this elimination of £ the dynamic equation (jl.lOj) is eliminated 
also on the foundation that the remaining Eulerian system of dynamic equations 
(|1.6|) - (|1.8|) for five dependent dynamic variables v, p, S is closed. Conventionally 
most researchers use RID, and hence the Eulerian system of dynamic equations. 
Another way of simplification is a use of the relabeling group for integration of 
the complete system of dynamic equations (jl.6|) - (jl.8J) . (jl.lOJ) and application of the 
system (fTTTTHl. (EHTT1) . ipjhjl . 

In general, equivalency of the system (jl.l(J|) . (|2.21|) . (|2.22j) and the system (jl.fcij) - 
(jl.8|) . (jl.lOj) can be verified by a direct substitution of variables p, S, v, defined by 
the relations (j2~THl) - (l2~27H) . into the equations <fT7B j> — (fTT5 |> . Using equations (fnUjl . 
(J2.21J) . f)2.22|) . one obtains identities after subsequent calculations. But such com- 
putations do not display a connection between the integration and the invariance 
with respect to the relabeling group (J3.4|) . Besides, a meaning of new variables if, rj 
is not clear. 
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We shall use for our investigations the variational principle (J3.1J) . Note that 
for a long time a derivation of a variational principle for hydrodynamic equations 
(ll.6j) - (|1.8j) existed as an independent problem [HI El 13 El EH EH EJ ■ Existence of 
this problem was connected with a lack of understanding that the system of hydro- 
dynamic equations (jl.6|) - (jl.8j) is a curtailed system, and the full system of dynamic 
equations (jl.6J) - (jl.8|) . (jl.lOj) includes equations (jl.lUj) describing a motion of the 
fluid particles in the given velocity field. The variational principle can generate only 
the complete system of dynamic variables (but not its closed subsystem). Without 
understanding this one tried to form the Lagrangian for the system p.fjj) - (|1.8jl as a 
sum of some quantities taken with Lagrange multipliers, lhs of dynamic equations 
(J1.6J) - (jl.8J) and some other constraints were taken as such quantities. 

Now this problem has been solved (see review by Salmon 9]) on the basis of the 
Eulerian version of the variational principle (j3.1|) . where equations (jl.lOj) appear 
automatically and cannot be ignored. In our version of the variational principle 
we follow Salmon jU] with some modifications which underline a curtailed character 
of hydrodynamic equations (jl.6|) - (jl.8J) . because the understanding of the curtailed 
character of the system (jl.6|) - (jl.8|) removes the problem of derivation of the varia- 
tional principle for the hydrodynamic equations (jl.6j) - (ll.8|) . 

The starting point is the action (j3.1|) . We prefer to work with Eulerian descrip- 
tion, when Lagrangian coordinates (particle labeling) £ = {£ , £}, £ = {£ x , £ 2 , £ 3 } are 
considered to be dependent variables, and Eulerian coordinates x = {x°, x} = {t, x}, 
x = {x 1 ,^ 2 ,^ 3 } are considered to be independent variables. Here £ is a temporal 
Lagrangian coordinate which evolves along the particle trajectory in an arbitrary 
way. Now £ is a fictitious variable, but after integration of equations the variable 
£ ceases to be fictitious and turns into the variable ip, appearing in the integrated 
system (fTTUJ) . (|T2T|) . 

Further mainly space-time symmetric designations will be used, that simplifies 
considerably all computations. In the Eulerian description the action functional (j3.1|) 
is to be represented as an integral over independent variables x = {x°,x} = {t, x}. 
We use the Jacobian technique for such a transformation of the action (J3.1j) . 

Let us note that according to (J2.5|) the derivative (J3.2|) can be written in the 
form 



dx a dJ ( dJ 



-i 



va = ^T = xT- 7«— ' a = 1.2,3. (3.6) 

Then components of the 4-flux j = = {p, pv} can be written in the form 

(j2.16J) . provided the designation (j2.19|) 

is used. 

At such a form of the mass density p the four-flux j = {j 1 }, i — 0, 1, 2, 3 satisfies 
identically the continuity equation (j2.11|) which takes place in virtue of identities 
flUD , d2~Hl) . Besides, in virtue of identities (|2~7}l . (I2~K1) the Lin constraints lOTty are 
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fulfilled identically 

dJ 

?d£ a = m p Q ^)—dii a = 0, a = 1,2,3. (3.8) 

a SO,i 

Components j l are invariant with respect to the relabeling group ([3.4)1 . provided the 
function p ($) transforms according to (|3.5|) . 
One has 

p Q (£)dtd£ = Po(^)^7 — dtdx. = —dtdx. (3.9) 
^so,o m 

m ( dx a \ 2 m ( dJ \ 2 ( dJ \ 2 



2 V dt J 2 \d^J V^W ' (3 ' 10) 

and the variational problem with the action functional ([3.1)1 is written as a varia- 
tional problem with the action functional 



1 / d J \ 2 ( d J \ 2 



E}pdtdx, (3.11) 



where p is a fixed function of £ = {£ , £} and of £ ai = d^ a , a = 1, 2, 3, i = 0, 1, 2, 3, 
defined by ()3.7|1 . is the internal energy of the fluid which is supposed to be a fixed 
function of p and So(£) 

E = E{p,S (£)), (3.12) 

where p is defined by ([3.7)1 and «So(£) is some fixed function of £, describing initial 
distribution of the entropy over the fluid. 

The action ([3.11)1 is invariant with respect to subgroup Qs °f t ne relabeling group 
(13.4)1 . The subgroup Qs is determined in such a way that any surface So(£) =const 
is invariant with respect to Qs - In general, the subgroup Qs is determined by two 
arbitrary functions of 

The action (j3.11[l generates the six order system of dynamic equations, consisting 
of three second order equations for three dependent variables Invariance of the 
action ([3.11)1 with respect to the subgroup Qs allow one to integrate the system of 
dynamic equations. The order of the system becomes reduced, and two arbitrary 
integration functions appear. The order of the system reduces to five (but not to 
four), because the fictitious dependent variable £ ceases to be fictitious as a result 
of the integration. 

Unfortunately, the subgroup Qs depends on the form of the function So(£) and 
cannot be obtained in a general form. In the special case, when So(£) does not 
depend on £, the subgroup Qs coincides with the whole relabeling group Q, and the 
order of the integrated system reduces to four. 

In the general case it is convenient to introduce a new dependent variable 

S = S (Z). (3.13) 

Addition of the new variable increases the number of dynamic variables, but at the 
same time this addition makes the action to be invariant with respect to the whole 
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relabeling group Q. It allows one to integrate the dynamic system and to reduce 
the number of dynamic variables. According to ()3.8j) the variable S satisfies the 
dynamic equation (J1.8)) 

fdiS = 0. (3.14) 

In virtue of designations (J2.16)) and identities ()2.7j) . (J2.8j) the equations (J3.14)) . (J3.8j) 
are fulfilled identically. Hence, they can be added to the action functional (|3.11|) as 
side constraints without a change of the variational problem. Adding ()3.14|) to the 
Lagrangian of the action 1)3.11)) by means of a Lagrange multiplier rj, one obtains 

* B - 51 = / {§ G£) ' (S " - pE + "MrM dti * (315) 

where E = E(p,S). The action 1)3.15)1 is invariant with respect to the relabeling 
group Q which is determined by three arbitrary functions of £. 

To obtain the dynamic equations, it is convenient to introduce new dependent 
variables j l , defined by (j2.16j) . Let us introduce the new variables j l by means of 
designations (j2.16j) taken with the Lagrange multipliers pi, i = 0,1, 2, 3. Considera- 
tion of ()2.16|) as side constraints does not change the variational problem, because 
conditions ([2.16)) are always compatible with dynamic equations, generated by the 
action ()3.15)) . Then the action ([3.15)1 takes the form 

A E [p,j,Z,P,V,S} = J^t-pE-pkfr-mpoiQ^^+ru'dkS^dtdx (3.16) 



For obtaining dynamic equations, the variables p,j,£,p,r),S are to be varied. 
Let us eliminate the variables pi from the action (j3.16[) . Dynamic equations arising 
as a result of a variation with respect to £ a have the form 

= t a p = -md k (p (g) ^7~H7 — Pi J + m -^T^^—Pk = °) a = 1, 2, 3 

(3.17) 

where C a are linear operators acting on variables p = {pi}, i = 0,1,2,3. These 
equations can be integrated in the form 

Pi = 9°(Qd£ + 9 a (mZ«, 2 = 0,1, 2, 3, (3.18) 

where £ is some new variable (temporal Lagrangian coordinate), g a {£), a = 1,2,3 
are arbitrary functions of the label g (£ ) is an arbitrary function of £ . The 
relations ()3.18)) satisfy equations ()3.17)) identically. Indeed, substituting ()3.18)1 into 
(J3.17)) and using identities ()2.6)) . ()2.7)) . we obtain 



- md k \ p (£) 



+ ^V(g = 0, a = 1,2, 3, 

(3.19) 
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Differentiating braces and using identities (|2.8j) . (|2.7|) . one concludes that ()3.19|) is 
an identity. 

Setting for simplicity 

dk<P = 9°(Qd k £ , A: = 0,1,2,3 (3.20) 

we obtain 

Pk = d kV + g a {£)d£ a , k = 0, 1, 2, 3 (3.21) 

Note that integration of equations (J3.17j) by means of the Jacobian technique 
and appearance of arbitrary functions g a (£) is a result of invariance of the action 
with respect to the relabeling group (|3.4|) . 

Substituting ()3.21|) in (|3.16|) . one can eliminate variables p i: i = 0,1,2,3 from 
the functional (j3.16|) . The term g a ($,)d k t; a dJ/dt; 0ik vanishes, the term d k (pdJ/d£ ok 
makes no contribution to the dynamic equations. The action functional takes the 
form 

Ag[p,j,t,V,S\ = J ^-pE-j k (d kV + g a ($)d k Z a -rid k S)^dtdx (3.22) 

where g a (£) are considered to be fixed functions of £ which are determined from 
initial conditions. Varying the action ()3.22|) with respect to <p, rj, S, j, p, we 
obtain dynamic equations 

6<p: d k j k = 0, (3.23) 
6£ a : n^j k d k ^ = 0, a = 1,2, 3, (3.24) 
where Q af3 is defined by the relations 

" = — , at, p = 1,2,3 (3.25) 

5r] : j k d k S = 0, (3.26) 
BE 

6S : j k d kV = -p—, (3.27) 

Sy. v=3/p = V<p + g a (t)VZ a -riVS, (3.28) 

6p : ~2?~ - doip ~ + ^ = °' (3 - 29) 

Deriving relations (|3.24|) . (|3.27j) . the continuity equation (|3.23|) is used. It is easy to 
see that ()3.24|) is equivalent to the Lin constraints (|1.1U|) . provided 

det || Q a/3 ||^ (3.30) 

If condition (J3.30J) obtains, and equations (|1.1U|) are satisfied, the equations ()3.26|) 
and ()3.23|) can be integrated in the form of ()2.18|) and ()2.19j) respectively. Equations 
(l3~2Tf) and (l3~2"51) are equivalent to flPD) and (l2~2"Uj) . Finally, eliminating d £ a and 
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8qS from (J3.29)l by means of (|1.10|) and (|3.26jl . we obtain the equation (J2.21)l and, 
hence, the system of dynamic equations (jl.lUj) . ()2.21|) . ()2.22|) . where designations 
(pH8j) - (f2~2(H ) are used. 

The curtailed system p.6jl - p.8j) can be obtained from equations (J3.23j) - (J3.29j) 
as follows. Equations (|S2SI), flUZED coincide with (fTfij) . IfTBjl . For deriving $TJ$ let 
us note that the vorticity uj = V x v and vxw are obtained from ()3.28|) in the 
form 

u = V x v = -fi^Vf^ x V£ a - Vr) x VS (3.31) 

v x w = H a,? V^(vV)( a + VS(vV)r/ - Vr/(vV),S (3.32) 

Let us form a difference between the time derivative of 1)3.28)1 and the gradient 
of ()3.29j) . Eliminating £l a Pdo£ a , d S and d rj by means of equations ()3.24jl . 1)3.26)1 . 
(|3.27)1 . one obtains 

- fi Q/3 V^(vV)£ a + V??(vV)S - VS(vV)r/ = (3.33) 
Using ()3.31|) . ()3.32)1 . the expression (|3.33)1 reduces to 

<V + V — + -V(p 2 — )-vx (Vx v) = (3.34) 
2 p ap 

In virtue of the identity 

v 2 

v x (V x v) = Vy - (vV)v (3.35) 

the last equation is equivalent to 1)1.7)1 . Note that at derivation of the curtailed 
system (11.6)1 - p. 8)1 the condition (J3.3()jl is not used, and the system (ll.fi)) - p. 8)1 is 
valid in any case. 

In the general case ()3.30)1 . differentiating equations ()3.28j) . ()3.29jl and eliminating 
the variables <p, 77, we obtain the curtailed system p. 6)1 - p. 8)1 . whereas the system 
(frmj) . (JOII), (l2~^2ll follows from the system (l3~25j) - (l3~2Hj) directly (i.e. without 
differentiating). It means that the system p.llljl . 1)2.21)1 . 1)2.22)1 is an integrated 
system, whereas the curtailed system p. 6)1 - p. 8)1 is not, although formally they 
have the same order. 

The action of the form ()3.22j) . or close to this form was obtained by some au- 
thors [Ulini, but the quantities g a , a = 1,2,3 are always considered as additional 
dependent variables (but not as indefinite functions of £ which can be expressed via 
initial conditions). The action was not considered as a functional of fixed indefinite 
functions g a (£). 

Thus, five equations p. 10)) . ()2.21)1 . ()2.22)1 with S, p and v, defined respectively 
by 1)2.18)1 . ()2.19j) and (J2.20)) . constitute the fifth order system for five dependent 
variables f = {£ ,£}, 77. Equations ([TBI. (Q)l . (H~TUH . flZZIMZZZ) constitute the 
seventh order system for seven variables p, <p, 77, S. 
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If Q al3 = 0, it follows from (l3~23J) that 



g a (£) = — $ (|) (3.36) 

Then it follows from (g^BD and (jORjl 

v = j/p = V (<p + <t> VS, (3.37) 

In the case of isentropic flow (VS* = 0) the quantities f2 a/3 coincide with vorticity, 
and the fluid flow is irrotational, as it follows from ()3.30j) . In this well as 

at fulfillment of ()3.37j) the dynamic equations ()3.24j) are satisfied due to relations 
= an( j Tj n cons traints ()1.10|) do not follow from dynamic equations (|3.24|) . 
In this partial case an addition of the Lin constraints (jl.lOj) to the curtailed system 
Q1.6j) . ()1.8p is not necessary. 

The system of equations ()3.23j) - (|3.29jl . as well as the system (jl.lOj) . (|2.18j) - 
1)2.22)1 contain full information on the fluid flow in the infinite space V. The system 
of equations (jl.lOj) . ()2.18)) - ()2.22)) is a system of five partial differential equations 
for five dynamic variables (p, rj, and one needs to give initial data for them. But 
the initial values for variables <p, rj, £ can be given in the universal form, which is 
the same for all fluid flows. For instance, one can set 

|(0, x) =£ in (x) =x, p(0,x) = 0, r/(0,x)=0 (3.38) 

Then according to ()3.28j) the initial values of variables p, S, v have the form 

v(0,x) = v in (x) = g(£)=g(x), (3.39) 
p(0,x) = p in (£) = p in (x) , 5(0,x) = 5 in (O = 5 in (x) (3.40) 

where p in (x), Si n (x), v in (x) are given initial values which determine the fluid flow. 
Variables £ label the fluid particles, and a choice of the form of labeling is unessential. 

Let now the form of <p(0,x) and T](0, x) be changed and, we have <p(0,x) = 
<p in (x), T] (0,x) = r] in (x), where <p in (x), and i] in (x) are some given functions. This 
replacement can be compensated by the change of g (x) 

g (x) -> v in (x) - V(f in (x) + 7]- m (x) VS in (x) 

in such a way that the initial values of p, S, v remain the same. All this mean that 
dynamic equations (jl.lOjl . ()2.18|) - ()2.22)) contain complete (full) information on the 
fluid flow, because the choice of initial values for variables ip, r], £ is unimportant 
for calculation of the fluid flow in the infinite space V. 

This is valid also for the fluid flow in the space region x 3 > 0, as it is shown in 
[TUj . In this case the boundary conditions for v at the boundary a; 3 = are expressed 
via the arbitrary functions g (£), £ = x G {x|x 3 < 0}, whereas the initial condition 
are expressed via the arbitrary functions g(£), £ = x G {x|x 3 > 0}. Apparently, 
functions g (x) determine the initial and boundary conditions for the velocity v also 
in the case of the fluid flow in any finite volume V, although it is not yet proved. 
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Thus, the information, which is essential for the fluid flow determination, is 
described by functions p in (x), Sj n (£), g . This information is introduced in dynamic 
equations in the form of arbitrary functions. Unessential information concerning 
the methods of the fluid description, is given by initial and boundary conditions 
for variables (p, rj, The variables £, (p, rj are auxiliary variables which represent 
a method of the fluid flow description. The Clebsch potentials £ label the fluid 
particles. At the same time the variables £ describe displacement of the fluid along 
their trajectories, and this description does not depend on the method of labeling. 
Variables p, rj describe separation of the velocity field v into parts. The form of this 
separation is inessential. The auxiliary variables £, <p, r\ are described by partial 
differential equations, whereas the fluid flow in itself is described by finite relations 
containing arbitrary functions p in (£), S in (£), g (£) (at fixed variables ip, rj, £, given 
as functions of t,x). 

Situation is rather unexpected and unusual. One can obtain a result of calcula- 
tion of the fluid flow, but it refers to the method of description (variables £, ip, 77), 
and this method of description is determined by partial differential equations. It is 
not clear how to resolve and to use this situation effectively. 

4 Variational principle for incompressible fluid 

Inviscid incompressible fluid of constant density is a special case of ideal fluid, when 
one may set p = p =const. In this case the continuity equation ()3.23j) . containing 
time derivative of density p and determining time evolution of p, becomes 

Vv = (4.1) 

Equation ()4.1|) does not contain time derivatives at all. It is rather a constraint 
imposed on initial values of velocity v, than a dynamic equation, describing evolution 
of one of dynamic variables. Formally it means that the system of hydrodynamic 
equations ceases to be hyperbolic and becomes elliptic. This circumstance changes 
the statement of the fluid flow problem. 

Let us set p = p = const in the action (|3.!22j) and introduce new variables 

v=j/p > p = const ( 4 - 2 ) 

It is easy to verify that rj = ??(£), S = S (£), and the last term of (|3.22|) can be 
incorporated with the term j k g a {£)dk£ a - Thus, the action for the incompressible 
fluid looks as follows 

2 

^ E [v, & y>] = Po J{\ - W<p - g a (£)d Z a - g a (£)vVUdtdx, (4.3) 

where g a (£) are arbitrary fixed functions of 
Variation with respect to v, £, ip gives 

5v: v = V<p + 9 a (£)V£ a (4.4) 
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Po ^ = + W^) = 0, a = 1, 2, 3 (4.5) 

Po 1 ^ = Vv = (4.6) 

In the general case the condition (|3.30J) is satisfied, and the multiplier Q af3 in 
(14. 5 j) may be omitted. 

Substituting (|4.4|l into (|4.5jl and (|4.fjjl . one obtains 

n * (dotp + (Vcp + <f (£) Ve 7 ) V^) =0, a = 1, 2, 3 (4.7) 
W + <? a,/3 (0 V^Ve Q + <f(£) V 2 £ a = 0, g^ = JL. (4.8) 

The dynamic equation for (p does not contain temporal derivative. If Q al3 = 0, the 
fluid flow is irrotational and dynamic equations (J4.7)) are fulfilled independently of 
the Lin constraints, which have the form 

dotp + (Vip + <f (£) V£ 7 ) = 0, a = 1, 2, 3 (4.9) 

Lin constraints (|4.9|) are not dynamic equations in this case. 

Conventional hydrodynamic equations for the incompressible fluid 

Vp 

Vv = 0, d v + vVv = (4.10) 

Po 

are obtained from relations (|4.4j) - ([4.6)) . Differentiating (|4.4|) with respect to t, we 
obtain 

9 v = V[d cp + g a U)d U - n^d ^V£ a (4.11) 

where is defined by (l3~23|) . It follows from (pL"I)l 

vx(V x v) = ^(OV^(W)^. (4.12) 

In virtue of ()4.5|) the last term in rhs of (|4.11|) coincides with rhs of ()4.12j) . Then 
using the identity (|3.35jl one obtains 

9 v + (vV)v = V (d 0( p + g a {£)d Z a + ^v 2 ) (4.13) 

The equation ()4.13j) coincides with the second equation (|4.1()|) . provided we use 
designation 

- = - - ^v 2 - d oV - g a (S)d Q Z a , p = const (4.14) 
Po Po 2 

Here the pressure p is determined after solution of the system of hydrodynamic 
equations (|Q |l -(jO) |l . or (ETTUD . 

Let us stress that the conventional form (J4.1Uj) of hydrodynamic equations is 
obtained from the hydrodynamic equations ()4.4j) - (j4.6|) by means of differentiation. It 
means that the form of hydrodynamic equations ()4.4j) - (j4.6|) is a result of integration 
of hydrodynamic equations (|4.10|) together with the Lin constraints (jl.lOj) . It is 
interesting also that the system of equations ()4.4|l - ([4.b)) contains time derivatives 
only in dynamic equations (|4.5j) . 
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5 Two-dimensional flow of incompressible fluid 



We have mentioned in introduction that statement of the problem of stream flow is 
different for the irrotational and rotational cases. This difference appears only after 
integration. The statement of the problem is different not only for two-dimensional 
flow of incompressible fluid. It is different for any inviscid fluid, and the source of this 
difference lies in the dynamic equations ()3.24|) (or 1)4.50 ) which exclude the degrees 
of freedom connected with rotation in the case Q al3 = 0. The method used for 
derivation of (J1.2)) cannot be used in general case. So in the following we present the 
method in a form which can be expanded to any ideal fluid. The method is based 
on a use of the generalized stream function (mutual application of the Jacobian 
technique and description in terms of potentials). To avoid technical complexities 
we apply this method to the case of two-dimensional flow of incompressible fluid, 
when the generalized stream function has only one component. 

Equations (jl.l)) and (jl.2j) are different. The first equation is linear, whereas 
the second one is quasilinear. The first one can be solved rather easily, whereas 
the second one can be solved only by means of an iteration procedure. Difference 
in complexity of equations (jl.lj) and (jl.2|) is technical, whereas the difference in 
statement of the problem is conceptual. The idea of our investigation is simple. 
Introducing the generalized stream function £ 2 we solve equations ()4.5|) . ()4.6|) for 
any £ 2 - Then from ()4.4)) we obtain dynamic equation for determination of £ 2 and 
necessary boundary conditions for this equation. 

Although it is possible to deal with equations ()4.4|) - (J4.6)) for incompressible 
fluid, we prefer to consider dynamic equations for slightly compressible fluid, whose 
internal energy has the form 



The incompressible fluid appears in the limit e — > 0. We consider dynamic equations 
obtained from the action ()3.22|) with the entropy S = const and rj = 0. These 
equations have the form 




E 0: p = const, 



e < 1 



(5.1) 









(5.2) 
(5.3) 







where Q a/3 is defined by the relation ()3.25|) . 



J 



/i= 1,2, 



(5.4) 



P 




/i, a = 1, 2 



9{pE) 
dp 




(5.5) 
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Here and in what follows a summation is produced over repeating Greek indices (1, 2) 
and over repeating Latin indices (0 — 2). We use designations (x°, x 1 , x 2 ) = (t, x, y). 
Comma before index k denotes differentiation with respect to x k , for instance, £ a k = 

Let us set 

j dJ^Q k d.J d(x^ 2 ) k = Q ^ 

d (x° , x 1 , x 2 )'' ' d£ ok d (x° , x 1 , x 2 ) ' 

where £ 1; £ 2 are some functions of (t,x,y). Equations (|5.2)L (j5.3|) are satisfied by 
()5.6|) for any functions £i,£ 2 > because of the identities 

dJ dJ 
dk 7^ = °> ^* = °. a = 1 > 2 ( 5 - 7 ) 

Substituting (|5.1|) in (|5.5j) . we obtain 

w( 1 + ;)'fc) = (-6- fl * + » B(a -")' (5 ' 8) 

In the limit e — > equation ()5.8|) becomes 



(5.9) 



Po Po^O^x 2 ) 

In the case of stationary flow we can assume that £ 2 does not depend on the time 

t 

Z x = Zi(t,x,y), i 2 = i 2 {x,y) (5.10) 

and 

£ 2 ,o = 0. (5.11) 

Equation (|5.9|) takes the form 

^1,1-^1,2-^0 = 0. (5.12) 

Let us consider equation (|5.12|) as an equation for dependent variable ^ with a given 
function £ 2 = £ 2 (x,y). The first order partial differential equation is equivalent to 
the system of ordinary equations 

dx dy dx _ d^ 

S2,2 S2,l S2,2 Po 

It is easy to see that 

C 2 = i 2 {x,y) (5.14) 

is an integral of the first equation (|5.13jl . Resolving equation (|5.14j) with respect to 
y, one obtains 

y = F 2 (x,C 2 ), £ 2 (x,F 2 (x,C 2 )) = C 2 (5.15) 
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Substituting ()5.15|) in the second equation (|5.13|) . one can integrate it in the form 



fi = Po 



X' 

Ci (t,C 2 ) + J 



dx 

£ 2 , 2 (x,F 2 (x,C 2 )) 



(5.16) 



where C\ is an arbitrary function of arguments (t, C 2 ). Here x = —L is the inflow 
surface S- m , where the inflow conditions are given. The inflow surface is placed in 
the fluid, and all dynamic equations must be satisfied on (S in . 

Differentiating the second equation (|5.15jl with respect to C 2 , we obtain 



£, 2 (x,F 2 (x,C 2 ))F 2 , 2 (x,C 2 ) 



(5.17) 



Here and in what follows the symbol ';' denotes differentiation with respect to one 
of arguments of the function. Index after symbol ';' denotes the ordinal number of 
the argument. By means of (J5.17j) integral in (J5.16)) may be written in the form 



dx 







e 2 , 2 (x,F 2 (x,C 2 )) dC 2 
The relation f|5.16jl takes the form 

d • 



J F 2 (x, C 2 ) dx 



(5.18) 



£i = Po 



C l (t,C 2 ) + — J F 2 (x,C 2 )dx 



(5.19) 



Differentiating (|5.1fjjl with respect to t, we obtain 

~0Ci(*,C 2 



Ci,o = Po 



Of 



PoCi;i (*> Q 



(5.20) 



JC 2 =5 2 (x,y) 

According to (|5.1U|) and (|5.6|) we have 



J 



d (x°, x 2 ) 
and equations ()5.4|) take the form 



2,2) 



2 _ d(£i,£ 2 ) _ 

J q t_ n _i \ — S1,0S2,1 



d (x°, x 1 ) 



(5.21) 



u = v 



Po 



di<p + 



2 £l, 0^2,1 a . a/t-\t 



Po 

The compatibility condition of equations ()5.22|) and ()5.23|) has the form 
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(5.22) 
(5.23) 



O ^1,0^2,2 ^1,0^2,1 



Po 



Po 



^ = |£ (5.24) 



/3 



Substituting (|5.20j) and using designation (J3.25)) we obtain 

^,ii + ^ 2 2 = -Po^ 12 (^i^ 2 ), (5-25) 
where the stream function %p is defined by the relation 

i; = C (t,Q = fc 1;1 {t,£ 2 )d£ 2 (5.26) 

Relation (J5.26)) is chosen in such a way that equations ()5.22|) . ()5.23|) coincide with 
(jl.4)) . Let us resolve relation ()5.26|) with respect to £ 2 in the form £ 2 = Q(t,tp), 
where Q (t, iji) is the functional inverse of Co (t, £ 2 ) 

C (t,Q(t,il>)) = i/; (5.27) 

Substituting £ 2 in ()5.16|) . we obtain Clebsch potentials £ a and £ 2 expressed via tfj by 
means of relations 



Ci (t,C 2 ) + — I F 2 (x,C 2 )dx 



C 2 =Q(t,ip) 

(5.28) 

According to (f5~2"Tj) . (j5^Hjl and fiTZty we obtain 

j 1 = -Po*P,2, f = Po^,i ( 5 - 29 ) 

Dynamic equation ()5.25)) for the stream function if) is the second order partial 
differential equation of elliptic type. It contains indefinite functions f2 12 (£i, £ 2 ), 
Co {t, C2) which are determined from the inflow conditions. The time t appears in 
dynamic equation only as a parameter describing evolution of the Clebsch potential 

The dynamic equation ()5.25|) is applicable for arbitrary two-dimensional flow of 
ideal incompressible fluid. 

Remark. Equation 1)5.25)1 can be obtained by direct integration of Euler equations 
for two-dimensional flow of incompressible fluid. Indeed, eliminating pressure from 
equations (J4.10)) for two-dimensional flow, we obtain 

u ,i + v )2 = 0, d 2 (w + m«,i + uu 2 ) = di (v fi + uv t i + vv )2 ) (5.30) 

where u and v are velocity components along rc-axis and y-axis respectively. Intro- 
ducing the stream function ijj 

«=-'0 J 2) v = ^,i ( 5 - 31 ) 

we satisfy the first equation ()5.30)) identically. The second equation 1)5.30)) takes the 
form 

d u - = °> u = tn + ^,22 (5-32) 

d(x,y) 

In the stationary case, the result of integration is considered in introduction. 
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In the case of arbitrary flow the equation (J5.32J1 is reduced to the form 



— = d uj + udiuj + v d 2 uj = 
at 



(5.33) 



Solution of ()5.33|) can be presented in the form 



LO = lf) >n + i>- 
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n(ei,e 2 ) 



(5.34) 



where £ x and £ 2 are two independent solutions of equation (|5.33|) . i.e. £ x and £ 2 
are Clebsch potentials, satisfying equations (jl.lOj) . Thus, the problem of the Euler 
system integration appears to be connected with consideration of the Lin constraints. 
The statement that we can investigate the fluid motion, considering only the Euler 
system of dynamic equations and ignoring completely Lin constraints, is valid only 
in the case Q = 0. In the general case it is only an illusion. In reality the Lin 
constraints are hidden inside the Euler system and appear after integration. 

In the simplest case, when all fluid particles cross boundary, and there are no 
closed stream lines, the function Q (^1,^2) ^ s determined from the inflow conditions 
on the boundary of the considered flow. If the flow contains fluid particles which 
move all the time inside the considered volume and do not cross its boundary, the 
function Q (£1, 62) i s t° be given for these particles also, although such an assignment 
is impossible via boundary conditions. Here we shall consider only the simplest case 
of rotational flow, when there are no closed stream lines. The case with closed 
stream lines of rotational flow inside the irrotational flow is considered in the paper 
|2Uj . In this paper the function Q (£ 2 ) = / (VO is given on the closed stream lines 
without a reference to boundary conditions. 

In the considered case the dynamic equation (J5.34|) appears to contain informa- 
tion on boundary conditions. It means that, studying rotational nonstationary flows, 
a separate investigation of dynamic equations and boundary conditions becomes im- 
possible. Appearance of several indefinite functions of if) in dynamic equations of 
type (jl.2j) for stationary rotational flow is well known (see, for instance, survey (3]). 
Unfortunately, an investigation of the origin of these functions is sometimes absent. 

Obtaining the basic equation (|5.25|) . we used properties of Clebsch potentials, 
described by equations (|5.4|) . These relations contain arbitrary functions g a (£), 
which describe influence of boundary conditions on the form of dynamic equations. 
Indefinite functions g a (£) describe general properties of ideal fluid, which could be 
applied not only for investigation of incompressible fluid. Information contained in 
three functions g a (£) of three arguments ^1,^25^3 i s m uch more, than information 
contained in one function Q (£ 1; £ 2 ) of two arguments. Maybe, not all this informa- 
tion is essential, but some part of this information is essential for sure. It is the 
reason why we investigate dynamic equations ()4.4|) - ()4.6|) . but not Euler system 
(J5.30|) . which contain this information only in implicit form. 

The last shorter derivation (J5.3U)) - (|5.34j) of equation (J5.34j) ignores general 
relations (|5.3|) . (|5.4|) . but nevertheless finally it leads to (|5.34|) . i.e. to a description 
in terms of Clebsch potentials. 
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6 Inflow conditions for stationary flow 



Let us suppose that the stationary inflow conditions are given at the plane «Si n : 
x = —L =const. (We consider the case, when there are no closed stream lines inside 
the considered volume.) The inflow surface Si n lies in the fluid completely, and all 
dynamic equations are satisfied on Si n . In the stationary case the stream function 
ip does not depend on time t. One may set in ()5.26|) Co (t,£ 2 ) = Co (£2)- On the 
inflow surface we choose the function £ 2 m the form 

&(-L,y) = y (6.1) 

According to (|5.21|) (J5.20J) the inflow conditions have the form 

u(-L,y) = -C 1 . >1 (-L,y) = V x (y) (6.2) 

^(-L,y) = -[di(^ (t,y)Ui(x,y))} =W(y) (6.3) 

OX Po ' /lx=-L 

where u and v are velocity components along rc-axis and y-axis respectively. V x (y) 
and W (y) are supposed to be known functions determining inflow conditions of the 
fluid. We obtain from (Q and (lo~2T)|) 

t 

Ci (t, y) = - fv x (y) dt = S (y) t, C (y) = - f V x (y) dy (6.4) 



The lhs of (I5.24|) at x = —L is equal to 

(-L, y)-^- y) = |-K (y) -W(y), (6.5) 
ay ox oy 

whereas rhs of ()5.24j) at x — —L is equal to p o r2 12 (^ 1 (t, — L, y) , y). Then using 
representation (j5.16j) . we obtain from (|5.24jl . ()5.28j) at x = —L 

p n 12 (p C 1 (t, y),y) = |-K (y) -W(y) = U (y) (6.6) 

where U (y) is known function of variables y. As far as rhs of (|6.6J) does not depend 
on t, lhs of ()6.6p does not depend on Ci, and Q 12 depends only on ^ 2 

^ 12 (^i, Q = ^ 12 (0, £ 2 ) = -U (£ 2 ) = - {-tL-Vx (Q ~ W (Q) (6.7) 

PO PO \^2 / 

Thus, inflow conditions (|6.2jl . (|6.3j) admit one to determine indefinite functions 
Q 12 and G\. It is interesting, that inflow conditions are given at the inflow surface 
S m : x = —L by the velocity component u and by derivative dv/dx (normal to «S in ) 
of the component v (but not by the component v itself). Information, giving by 
the velocity component v on the inflow surface <Sj n , is insufficient for determination 
of indefinite functions Q 12 and C\. The derivative dv/dx together with du/dy 
determine vorticity. The derivative du/dy can be calculated, if the velocity u is 
given on the inflow surface, whereas dv/dx cannot. To determine vorticity, we need 
to give dv/dx on the inflow surface, whereas the tangent component v itself appears 
to be unnecessary for the vorticity determination. 
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7 Problem of stationary flow over obstacle 



Let us consider a flow of incompressible fluid inside the space region V bounded by 
two walls S\ : y — —R =const, S% : y = R, by inflow surface S- m : x = —L =const, 
by outflow surface «S out : x = L, and by the surface iS b of some obstacle, placed 
in vicinity of the coordinate origin. The shape of the obstacle is determined by the 
relation 

S ob : S ob {x,y) = (7.1) 
The boundary conditions are formulated as follows. 



<Si : ip (x, —R) — if) i — const, S 2 : ip (x, R) = ip 2 = const 
S oh : ip (x, R) = ip oh = const 

Constants ■0 1 and ip 2 are connected by the relation 



(7.2) 



R 

V>2-V>i = -/ V x (y)dy (7.3) 

-R 

Conditions ()7.2|) mean that the normal velocity component vanishes on all walls. 
Value of ij) oh describes a circulation around the obstacle, because ip a b ~ ^1 describes 
the mean value of flux between S\ and iS b- For instance, if ip ob = if) 1 , the circulation 
around the obstacle is such, that the fluid flows only between *S 2 and <S b 
Boundary condition on the inflow surface iS in has the form 



S m : i>{-L,y)=i> 1 - J V x (y)dy (7.4) 

-R 

where the velocity component V x (y) along the x-axis is presumed to be given. The 
boundary condition on the outflow surface S out has the form 

Sout : iff {L, y) = ^out (y) , V'out = V'out ( R ) = ^2 (7-5) 

where ip ont (y) is some given function. It describes the velocity distribution on the 
outflow surface. To understand, why such a distribution should be given, and why it 
is not determined by the inflow condition, let us imagine that outside the volume V 
in vicinity of S ont there are some obstacles. Distribution of velocity on <S ou t depends 
on size, shape and position of these obstacles. Description of these obstacles is 
carried out by ^ out (y). 

Let us consider a flow over the obstacle and rewrite dynamic equation ()5.25|) in 
the form 

Cip = n(ij) (7.6) 
£V = V 2 ^ = ^,11 + ^,22 (7-7) 

n (vo = - Po ^ 12 (Q = -u (£ 2 ) = -u (Q m (7.8) 
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where the function U is defined as a function of argument y by means of the relation 

dnsD 

U(y) = ^V x (y)-W(y) 
Function £ 2 = Q (ip) is defined as a function of argument ip by means of the relation 

^=■01- / Kb (j/) rfj/ (7.9) 



-R 

In particular, if (y) = V =const, £ 2 = Q (VO — ijPi ~ 4>) /Vo- 

Dynamic equation ()7.6|) is a nonlinear equation, in general. Its solution is diffi- 
cult. The goal of this paper is to provide a statement of the stream flow problem, 
because even this problem is not solved for irrotational flows. So far as the solution 
of equation ()7.6|) is concerned, it is a special mathematical problem which can be 
solved by different methods. 

If vorticity Q is small, dynamic equation ()7.6|) can be solved by means of an 
iteration procedure. To obtain the zeroth order approximation V(o)> we consider 
linear equation 

£V( ) = V 2 ^ (0) = (7.10) 
with boundary conditions (|7.2|) - (|7.5|) . which is written for brevity in the form 

V'(O)] E = V>2 (*> X ) > XGE I 7 - 11 ) 

where E is the boundary of the volume filled by fluid. 

Let G ( x l x ') be the Green function for equation (J7.1Uj) . It is defined by the 
relations 

V 2 G (x|x) = -4tt5 (x - x') , G (x|x) ^ = (7. 12) 

Solution of equation (|7.1U|) is given by the formula 



1 r dG ( xlx') 

^(0) (*) = - T J -^r^^ OO ^ (7-13) 

s 

where d/dn' means derivative with respect to x' in the direction normal to E. Solu- 
tion ()7.13|) is an irrotational approximation to the solution of ()7.6|) for the rotational 
flow. 

The first approximation is determined by the relation 

4, {1) (x) = ^(o) (x) - ^ / G (x|x') n (V» (x')) dV (7.14) 

v 

where integration is produced over the volume V filled with the fluid. 

iP (n) (x) = ^ (0) (x) - i- I G (x|x') O (x')) dV', n = 2, 3, ... (7.15) 

v 
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If the iteration process converges, we obtain a solution of the problem of rotational 
stationary flow around an obstacle. One can hope that the iteration process appears 
to be effective at least for small vorticity Q. 



8 Examples of two-dimensional stationary flow. 

Let us consider a stationary flow, bounded by two parallel plane walls Si : y — 
0, £2 : y = b, with inflow and outflow surfaces iS; n : x = 0, S ont : x = a. 
To obtain linear equation ()7.6j) which has a simple exact solution, we choose the 
following boundary condition 

ip(x,0) = 0, if>(x 1 b) = 0, (8.1) 

7iy Try 
ip(0,y) = Vosin— , ip (a, y) = V sin — , V Q = const (8.2) 



ip u (0,y) = fcVosin^p k= y Q+ (^)> ft = const (8.3) 

In this case the vorticity Q (i/j) in the rhs of ()7.6|) on the boundary x = is determined 
by the first condition (Q and (jFOj) 

JLu(0,y)-±v(0,y) = ^ 11 (0,l/) + ^ 22 (0,l/)=(/c 2 -(^) 2 )K ) sin^ 

= W(0,y) 

It means that the function Q (-0) in the rhs of (J7.6|) has the form Q (ip) = flip, Q =const, 
and dynamic equation takes the form of linear equation 

V 2 ^ = Qip, Vl = const (8.4) 

Its solution, satisfying the boundary conditions ()8.1|) . (|8.2jl has the form 



Tr sinh (kx) + sinh (A; (a — x)) . ny . L /tt\ 2 

This solution describes a set of flows with vanishing total flux through the inflow 
surface and different vorticities Q. Flows of this set are distinguished by the power 
of circulation which is described by the parameter Q. 

Let us consider two-dimensional stationary slightly rotational flow between two 
concentric circular cylinder of radius a and b (a < b). Let us choose the coordinate 
system with the z-axis directed along the cylinder axis and coordinate origin on 
the cylinder axis. The x-axis is directed along the mean velocity Vo of the flow. 
Let us introduce polar coordinates (r, ip) in the (x, y)-plane and set the boundary 
conditions for the stream function as follows. 

*j;(a,<p) = 0, VKM = -^o(l-^Wn^ (8.6) 



1 d I dip (r, (p) ~ 
r dr \ dr 



r=b 



^>=«W<M) = {t> I |*|^(8.7) 
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where Ui (ip) is some given function of ip. It is supposed that the flow is rotational 
only in the part of fluid which flows near the cylinder. Let — ipo and i/j 

A = V (l- Na (8.8) 

be values of the stream function on the boundary between the rotational and irrota- 
tional parts of the flow. Na is the transverse size of the rotational region (N 3> 1). 
After statement of the problem we set b — > oo. The irrotational flow of the zeroth 
approximation around the cylinder is described by the stream function ip^ which 
is given by the expression [I] 



V>(o) ( r >V?) = -V il- — J rsin^, 



where r, ip are polar coordinates on the (x, ?/)-plane 

x = r cos ip, y — r sin ip 



r > a 



.9) 



(8.10) 



The irrotational flow (J8.9)) is used at formulation of boundary conditions on the 
external cylinder of radius b. 

According to (|7.14j) and ()8.6|) - ()8.9|) we obtain for the first approximation 



^(i) (r, <p) = ip (0) (r, ip) - 



47T 



G (r, ip\r', ip') uj (-0(0) (r', <p'fj r'dr'dip' 



\-4>(o)( r 'y)\<^o 



(8.11) 

where u (ip) is defined by the relation (|8.7J) . and the Green function G (r, y>|r', y/) is 
taken for the case b = oo. It has the form 



G (r, <p\r', ip') = In 



a 2 + rL ^ 2rr' cos (ip — ip') 

r 2 + r' 2 — 2rr' cos (ip — ip') 



and integration is produced over the region, where V(o) i r i v) < =const. Ac- 
cording to ()8.8|) . ()8.9|) this region is defined by the relation 



Vo l r sin <p 



< V Na, 



r > a, b = oo 



where AT > 1. If 



A" 



r<i? = a — + \ — +1 



N 



(8.13) 



(8.14) 



the inequality ()8.13|) is satisfied for any angle ip, and for r' < R one should integrate 
in 1)8.11)1 over all values of angle ip' . For r' > R it is useful to introduce the variable 
?7 by means of relation 



ip' = $ (r, 77) = arcsin 



Vn 1 



(8.15) 
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and integrate over rj in the region, where r' > R. 
We obtain 

V'(i) = "0(0) + h + h + h (8.16) 

where 



R 2vr 



h = ~ jr'dr 1 ' J G{r^\r'^')u(i, m (t~V)) V (8.17) 

a 

h=— r dr / * (8.18) 

r 1 7 aj /^^^(r^ir'^Cr', 77)) a; (-77)^(1-^)^ 

3 = iw r y / 2 2 — — ( 8 - 19 ) 

% 7r-$(r'^ ) ^ 2 (l-^) r' 2 -77 2 

where -0 O = VoNa. Although calculations of integrals ()8.17|) - ()8.19|) is not simple, 
they give the first approximation of the rotational flow around the cylinder. 



9 Conclusion 

Theory of rotational flows is more complicated, than that of irrotational flows. There 
are two reasons for this. 

1. Dynamic equations for rotational flows appear to be nonlinear almost always. 
They cannot be solved exactly. 

2. The problem of rotational stream flow is not stated yet properly. 

The first problem is pure mathematical. Existence of this problem is quite clear 
for all researchers. The second problem is physical. We believe that hardly anybody 
of researchers guesses on existence of this problem. The first problem is very difficult. 
The second problem is not so difficult, but it is much more important that the first 
problem. Why? 

First, because existence of the second problem in itself seems to be problematic. 
Second, the irrotational flow is a special degenerate case of the rotational flow, when 
vorticity vanishes. Transition from the general case to the degenerate one is rather 
simple. The reverse transition from the degenerate case to the general one is much 
more difficult, because the general case contains additional information which is 
absent in the degenerate case. If we could not evaluate correctly this additional 
information, we do not enable to construct a perfect theory. Practically, we ignore 
dynamic equations (jl.lUJI and boundary conditions for determination of vorticity. 
As a result we obtain imperfect theory, where some degrees of freedom are ignored. 
In some cases these hidden degrees of freedom manifest themselves suddenly, and 
we try to understand what could be meant by this effect. At any rate a serious 
mathematical investigation of such difficult hydrodynamical problems as turbulence 
is impossible, if the underlying rotational flow theory is incomplete. 
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